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Abstract 



Context. The coming Planck and Herschel missions will survey the sky at unprecedented angular scales and sensitivities. 
Simulations are needed for better interpretating the results of the surveys and for testing new methods of, e.g., source 
extraction and component separation. 

Aims. We present new simulations of the infrared and submillimeter cosmic background, including the correlation 
between infrared galaxies. The simulations were used to quantify the source-detection thresholds for Herschel/SPIRE 
and Planck/HFI, as well as to study the detectability of the cosmic infrared background correlated fluctuations. 
Methods. The simulations are based on an empirical model of IR galaxy evolution. For these correlations, we only 
included the linear clustering, assuming that infrared galaxies are biased tracers of the dark-matter fluctuation density 
field. 

Results. We used the simulations with different bias parameters to predict the confusion noise for Herschel/SPIRE 
and Planck/HFI and the completeness levels. We also discuss the detectability of the linear clustering in Planck/HFI 
power spectra, including the foreground and backgrounds components. 

Conclusions. Simulated maps and catalogs are publicly available online at 
http: / /www. ias.u-psud.fr/irgalaxies/simulations.php 

Key words, infrared: galaxies -galaxies: evolution - (cosmology:) large-scale structure of universe 



1. INTRODUCTION 

The cosmic infrared background (CIB) (A > 8/ito) is the 
reHc emission of the formation and evolution of galaxies. 
The first ob servational e v idence of this background was 
reported by iPuget et al] (119961) and t h en co nfirmed by 
iHauser et all (|1998D and iFixsen et aP (|l998l ). The dis- 
covery of a surprisingly high amount of energy in the 
CIB has shown the importance of studying its sources 
to understand how the bulk of stars was formed in the 
Universe. Deep cosmol ogical surveys have been carried 
out t hanks to ISO (see iGenzel fc Cesarskvl l2Q00l : lElbaj . 
20051 for reviews) mainly at 15 with ISOCAM (e.g. 



Elbaz et all. 120021): at 90 and 170 fim with ISOPHOT 
(e.g. iDole et al l. I2OOII): to S PITZER at 24. 70, and 



160 /Im (e.g. iPapovich et all . 12004 iDole et al.1 . I2Q04D 
and to groun d-based instruments such a s SCUBA (e.g . 
Holland et all . [1998) and MAMBO (e.g. iBertoldi etall . 



200dh at 850 and 1300 ^m, respectively. These surveys 



have allowed f or a better understa nding of the CIB and 
its sources (see iLagache et aT I, l20n,l for a general review). 
Some of the results include: the energy of the CIB is 
dominated by starbursts although AGN (active galactic 
nucleus) contribute too, and the dominant contributors to 
the energy output are the LIRGs (luminous IR galaxies) at 
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z ^ 1 and ULIRGs (ultra luminous IR galaxies) at z ^ 2—3. 

Determination of the CIB by the COBE satellite 
has been hindered by the accuracy of subtracting the 
foreground b y only providi n g just upper Hmits at 12, 25, 
and 60 (iHauser et al.. Il998l). lower limit has been 
derived at 24 fj,m by Papovich et al.l (|2004f ) as well as the 
contributi on o f 24 fim galaxies to the background at 70 and 
160 /zm (Dole et al., 2006). The contribution of the galaxies 
down to 60 fiJy at 24 /im is at least 79% of the 24 fim 
background, and 80% of the 70 and 160 background. 
For longer wavelengths, recent studies have investigated 
the contribution of populations selected in the near-IR to 
the far- infrared background (FIRB, A > 200 iim): 3.6 urn 
select ed sources to the 850 /xm background (jWang et al.l . 
I2006h and 8 irm and 24 /im s elected sou r ces to the 850 
jim and 450 /im backgrounds ijDye et al.l . l2006f ). Similar 
studies with Planck and Herschel will provide even more 
evidence of the nature of the FIRB sources. 



Studying correlations in the spatial distribution of IR 
galaxies as a function of redshift is an essential observation 
(parallel to the studies of individual high- redshift, infrared, 
luminous galaxies), to understand the underlying scenario 
and physics of galaxy formation and evolution. A first study 
has been done using the 850 /um galaxies ijBlain et al.l . 
|20Q4( ). Although the number of sources is quite small, they 
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find evidence that submillimiter galaxies are linked to the 
formation of massive galaxies in dense environments des- 
tined to become rich clusters. This has now been directly 
supported by the detection of the clustering o f high-redshift 
24 /im selected ULIRGs and Hy perLIRGs (|Farrah et all 
I2OO6I : iMagliocchetti et all l2007t ). Studying correlations 
with individual IR galaxies is very hard due to either high 
confusion noises, instrumental noises, or small fields of 
observation. It has been shown that the IR-background 
anisotropics could provide information on the correlation 
between the sou rces of the GIB and dark matter for large- 
scale structures (jKnox et al.l . I2OOII : iHaiman fc Knoxl . l2000l 
hereafter HK) and on the large-scale structure evolution. 
First studies at long wavelengths have only detected the 
shot-n oise component of the fi uctuatiqns: Lagache fc Pugetl 
II2000I1 at 170 MTO, [Ma tsuhar a et al l II2OOOD at 90 and 170 
um. iMiville-Deschi^s et al.l ll200lt) at 60 an d 100 urn. 
iLagache et all l|2007t l and iGrossan fc SmootI l|2007l l re- 



port first detections of the correl ated componen t using 
Spitzer/MIPS data at 160 fim. iLagache etall (120071 ) 
measured a linear bias 6 ~ 1.7. 

Future observations by Herschel and Planck will allow 
us to probe the clustering of IR and submm galaxies. 
Nevertheless these experiments will be limited, the confu- 
sion and instrumental noises will hinder detections of faint 
individual galaxies. Clustering thus has to be ana l ysed i n 
the background fiuctuations (e.g. iNegrello et al.l . l2007l l. 
The need for a prior understanding of what could be done 
by these experiments has motivated us to develop a set of 
reaHstic simulations of the IR and sub-mm sky. 

In Sect. 2 we present the model on which are based 
our simulations. In Sect. 3 we discuss how the simulations 
are done and present a set of simulated sky maps and 
their corresponding catalogs. Different catalogs are created 
for 3 different levels of correlation between the IR galaxy 
emissivity and the dark-matter fiuctuation density field 
(strong, medium, and no correlation). For each of these 
catalogs, we can create maps of the sky at any given IR 
wavelength and simulate how different instruments will 
see them. We focus in this paper on Planck/HFI and 
Herschel/SPIRE. In Sect. 4 we use the simulated maps to 
give predictions for the confusion noise, the completeness, 
and the detection limits for each of the study cases, 
including the instrumental noise. In Sect. 5 we present the 
power spectra of the GIB anisotropics for Planck/HFI and 
discuss their detectability against the significant sources of 
contamination (shot noise, cirrus, and cosmic microwave 
background (CMB)). 

Throughout the paper the cosmological parameters were 
set to ft. = 0.71, f^A = 0.7, n„, = 0.27. For the dark-matter 
Hnear clustering we set the normahzation to ag = 0.8. 



2. THE MODEL 

2.1. Galaxies' empirical evolution model 

The model of IR galaxies us ed f or the simula- 
tions is fro r n ILagache et al.1 (120031 1. revisited in 
Lagache et al.1 120041 ) - hereafter the LDP model, see 
http://www.ias.u-psud.fr/irgalaxies/model.php, This 
model is a fiexible tool for planning surveys and developing 



analysis methods. The requirement was to build the 
simplest model of the luminosity function (LF) evolution 
with redshift, with the lowest number of parameters, but 
accounting for all statistical observational data between 5 
and 1 mm. These are the spectral energy distribution 
of the GIB and its fiuctuations, galaxy luminosity functions 
and their redshift evolution, as well as the existing source 
counts and redshift distributions. 

The luminosity function of IR galaxies was modelled by 
a bimodal star-formation process: one associated with the 
passive phase of galaxy evolution (normal galaxies) and 
one associated with the starburst phase, mostly triggered 
by merging and interactions (starburst galaxies) . Unlike for 
the starburst galaxies, the normal galaxy contribution to 
the luminosity function was considered mostly unchanged 
with redshift. The spectral energy distribution (SED) 
changes with the luminosity of the source but is assumed 
constant with redshift for both populations in this simple 
model. 

This model fits all the experimental data and has pre- 
dicted that LIRGs (10" < Lir < 10^^) dominate at 
z ~ 0.5 - 1.5 and that ULIRGs/HLIRGs {Lir > lO^^) 
dominate at z ~ 2 — 3 the energy distribution of the CIB. 
One example of the agreement between the model and the 
observations is shown in Fig. [2j 

2.1.1. Number counts and CIB fluctuations 

To illustrate the interest of studying the cosmic back- 
ground fiuctuations in the far-IR and submm domains, we 
use a simplistic approach for the number counts, following 
Lagache & Puget (200(|). The source number counts can be 
schematically represented by a power law: 



N{S > 5o) = iVo ( 

Do 



(1) 



where we set Sq to be the detection limit for the sources 
and iVo the number of sources with fiux larger than Sq . 

In a Euclidean Universe with uniform density of the 
sources a = 1.5. In the far-IR and submm, a steeper slope 
is observed with a = 2 — 3 in the regime where negative 
K-correction dominates. As an example , ISO observation s 
found a slope of a = 2.2 at 170 fj.m ijPole et all . l200lt ). 
Obviously, the number counts need to fiatten for low fiuxes 
to ensure that the GIB remains finite. For the rest of the 
discussion we will assume that a = for S < S*. The total 
intensity of the GIB composed by all the sources up to Smax 
is given by: 



IciB 



S—-—dS. 
ab 



For the Euclidean case the GIB intensity is dominated 
by sources near S* . 

Fluctuations from sources below the detection limit So 
are given by 



^" rIN 

dS 
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Figure 1. Emissivities computed using the LDP model at Figure 2. Comparison of the observed source counts (data 

(observed) 250 fim (continuous hne), 450 fim (long-dashed points) and model predictions (continuous lines). Upper left 

Hne), and 850 ^im (dotted-dashed line). The emissivity from 24 urn, Upper right 70 urn, Lower left 160 ^m, Lower right 

HK at 450 (short-dashed line) is shown for comparison 850 /im. 



Using ^ given by Eq. [T] we get 



2- a 



NoS^ 



For a > 2 CIB fluctuations are dominated by sources 
close to S* so that the same sources dominate both the 
FIRB and its fluctuations. Therefore by studying the fluc- 
tuations of the FIRB, we are also studying the sources that 
form the bulk of the contribution to the FIRB. We can 
check this conclusion with the number counts from the LDP 
model. Figure [3] shows that the same sources dominate the 
background and the fluctuations, but only for faint sources 
(for example S'sso < 50 mJy). Therefore, it is necessary 
to subtract bright sources prior to any fluctuation analysis 
since they would otherwise dominate the fluctuations. 

2.1.2. IR galaxy emissivity 

For the purpose of the model we need to compute the 
mean IR galaxy emissivity per unit of comoving volume 

[W/Mpc^/Hz/.i 



It is deflned as 

dN 



jd{v, z) = (1 

where L is the luminosity (in W/Hz/sr), ^^^^^^ ^, 



IS 

the comoving luminosity function (in Mpc~^), and v the 
observed frequency. We compute jd using the SEDs and 
luminosity function from the LDP model which assumes 
that the SED depends only on Lhoi. The resulting jd is 
different from wh at is used by former approaches (HK, 
iKnox et all . l200lf ). We can see the difference between the 
emissivity from our model and that of HK in Fig. [TJ The 
crude model used for the emissivities by HK gives much 
lower emissivities than ours. 



2.2. IR galaxy spatial distribution 

Any model trying to account for CIB fluctuations must 
describe the statistical properties of the spatial distribution 




-3 -2 

Log S(Jy) 

Figures. Contributions of the sources of flux S (in Jy) 
per Log interval of S to the background (dotted line, right 
y axis) and fluctuations (continuous line, left y axis) at 
850/Ltm. 



of the sources. The absence of a completely developed 
theoretical model for the distribution of the whole IR 
galaxy populations makes the empirical modelling of dif- 
ferent distributions for the sources of the CIB necessary in 
order to prepare future observations. We used an empirical 
description for the spatial distribution of these sources, 
which has been used to create the simulated sky maps. 

The LDP model did not address the spatial distribution 
problem due to the lack of constraints at the time it was 
built. This is still mostly the case at the t ime o f writing 
this work. The simulations by iDole et aP l|2003l ) did not 
implement any correlation between IR galaxies and used 
an uncorrelated random distribution. However, since future 
experiments such as Herschel and Planck will be able to 
detect large-scale IR galaxy correlations {i < 1000), a 
model addressing this problem has become necessary. 
Herschel, with its high angular resolution, is expected 
to also probe correlations between galaxies in the same 
dark-matter haloes but in this study this correlation has 
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Figure 4. Top: CIB Power spectrum with a bias 6=1 
at Herschel/SPIRE wavelengths 500 /im (continuous line), 
350 /im (dashed line), and 250 /im(dotted-dashed line). 
Bottom: CIB Power spectrum with a bias & = 1 at Planck/HFI 
CIB wavelengths 850 fim (continuous line), 1380 fim (dashed 
line), and 2097 /im(dotted-dashed line). 



not been considered for simplicity. We only consider the 
Hnear clustering i.e. IR galaxies as biased tracers of the 
dark matter haloes, with a Hnear relation between the 
dark-matter density-field fiuctuations and IR emissivity. 



We follow the prescription from lKnox et al.l (|2Q0lf ). The 
angular power spectrum that characterises the spatial dis- 
tribution of the fiuctuations of the CIB can be written as 

J '^^/{z)-jl{v,z)h\k,v,z)PM{k)\u=i/rG\z). (2) 

In the equation several components can be identified, 
starting with a geometrical one ^j^a?{z) (these terms take 
all the geometrical effects into account), followed by the 
galaxies emissivity jd{v, z) already described in Sect. 2.1.2, 
then the bias b(k,i>,z), and finally the power spectrum 
of dark-matter density fiuctuations today -P/\f (fc)|A;=;/r and 
the linear theory growth function G^(z). Finally I is the an- 
gular multipole, in the Limber approximation k = l/r, and 
r the proper motion distance. The way the power spectrum 
has been obtained is developed in the following subsections. 



2.2.1. Dark-matter power spectrum 

The power spectrum of the dark-matter distribution at z = 
can be written as 



PM(fc) oc kT^{k) 



(3) 



where T{k,t) is the transfer funct ion for a cold dark- 
matter universe (jBardeen et al.l . Il986l ) . The linear theory 
growth function G{z) writes as 



G'iz) 



g^n{z),nA{z)) 
g^ino,nAo)ii + zy 



(4) 



with 



g[n{z),nAiz)] = -n{z)x \n{z)'^^'^ - nA{z) 



And flA{z) 



n{z) 



no{i+zy 
Oo(i-H2)^+i-no 



, 2.2.2. Bias Model 



The bias of IR galaxies represents their level of correlation 
with the dark-matter density field. It can be expressed as 
a function of the spatial scale, the redshift, and the wave- 
length of observation. In this paper and due to lack of mea- 
surements for the bias, we consider a simplified constant 
bias b. 

Sjj{k,iy, z) ^ ^ 5p{k,iy, z) 
jd{k,iy,z) p{k,v,z) 

where jd is the emissivity of the IR galaxies per 
comoving unit volume, jd its mean level, and 5jd its 
fiuctuations. Similarly, p is the dark matter density, p 
its mean value, and 5p is the linear-theory dark-matter 
density-field fiuctuation. 

We have better knowledge of the bias for optical and 
radio galaxies than for IR galaxies. Several studies have 
been able to measure the bias for the optical sources. 
As an example, a high bias (6 ~ 3) has been found at 
z ^ 3 for the Lyman-Break Galaxies jS teidel et all. 119981 : 
iGiavalisco et al.1 . Il998l : lAdelberger"et al.. . ,199^ ))~^It has 
been found as well t hat the bias increase s with redshift 
both for the optica l ijMarinoni et al.l . l2006f ) and the radio 
(jBrand et al.l . l2003 l populations. The optical or radio bias 
could be misleading as a first guess for the bias of IR 
galaxies. IRAS has measured a low b ias of IR galaxies 
at z ~ (e.g. ISaunders et al.1 . Il992h . Such a low bias 
is expected since the starburst activity in the massive 
dark-matter haloes in the local universe is very small. But 
we expect a higher IR bias at higher z, du ring the epoch 
of for mation of galaxy clusters. Indeed, iLagache et al.l 
(|2007f ) report the first measurements of the bias, h ^ 1.7 
in the CIB fiuctuations at 160 pm using Spitzer data. The 
LDP model indicates that galaxies dominating the 160 
pm anisotropies are at z 1. This implies that infrared 
galaxies at high redshifts are biased tracers of mass, unHke 
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Figure 5. Redshift contributions to the angular power 
spectrum ^ at i = 1000 in i^K"^ for different wavelengths: 
250 fim (continuous line), 350 nm (dotted line), 550 fim 
(dashed line), 850 fim (dotted-dashed line), and 1380 /zm 
(long dashed line). 



in the local Universe. For an ext ensive review of the bias 
problem see iLahav fc Sutd (|2004j). 

The IR bias could have very complex functional 
dependences, namely with the spatial frequency k, the 
redshift z, and the radiation frequency v (for example 
if different populations of galaxies with different SEDs 
have different spatial distributions). However, for the 
simulations, simplified guesses for the bias were used, 
namely a constant bias of 1.5, 0.75, and 0. 

Figures |4] show the angular power spectrum C; for some 
Herschel and Planck wavelengths. The power spectra are 
shown for a constant bias 6=1. Since Ce oc 6^, the pre- 
dicted power spectrum scales as b^. 

2.3. Discussion and implications of the model 

The IR galaxy SED peaks near 80 /im. This combines 
with the Doppler shift and causes observations at different 
wavelengths to probe different redshifts. Figure [5] shows 
the contributions to the power spectrum at I — 1000 for 
different redshifts, normalized to unity. The contributions 
to the same £ come from higher redshift as wavelengths 
increase. The shorter wavelengths probe the lower redshifts 
because they are close to the maximum of the SED, while 
the longer wavelengths probe the higher redshifts due to 
the strong negative K-correction. 

Figures [7] and [6] show the redshift contributions to the 
intensity of the CIB and to its integrated rms fluctuations 
for Planck/HFI and Herschel/SPIRE, assuming sources 
with S > Sdet have been removed - Sdet corresponds to 
the source detection thresholds computed in Sect. 4.3. We 
see that the fluctuations and the FIRB are dominated 
by sources at the same redshift. Therefore, studying the 
fluctuations at different wavelengths will allow us to study 
the spatial distribution of the sources forming the FIRB at 
different redshifts. 
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Figure 6. Redshift contribution to the FIRB {top panels) 
and its fluctuations {middle panels). Also shown are the 
redshift distributions of the detected sources {bottom pan- 
els) for a typical large Herschel/SPIRE deep survey (see 
Sect 4.3). From top to bottom: 250/im, 350/iTO and 500/im. 



The amount of fluctuations that come from sources 
at redshifts lower than 0.25 for the Planck/HFI case at 
350 fim is noticeable from Fig. [7l This contrasts with 
the Herschel/SPIRE predictions where the bulk of the 
low-z sources contributing to the fluctuations in the 
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Figure 7. Redshift contribution to the FIRB {top panels) and its fluctuations {bottom panels) for a Planck simulation 
(dz=0.25) at 350/j,TO (top-left figure), 550/iTO (top-right figure), 850^to (bottom-left figure), 1380/im (bottom-right figure). 
The plots are for simulations with b = 1.5, which sets the detection limit (see Sect 4.3). 



Planck case are resolved. These individual detections 
with Herschel/SPIRE could allow their subtraction in the 
Planck maps. A similar approach could be used between 
the Herschel 500 fim and the Planck 550 /im channels, 
although it is more marginal. Using information on the 
fluctuations at shorter wavelengths to remove the low-z 
fluctuations from longer wavelength maps could be another 
approach to studying the fluctuations at high redshifts 
directly. 

A similar raodel h as been developed by HK and revisited 
bv lKnox et al.l l|200ll l. We compare the HK and our C'l pre- 
diction at 850 /im in Fig. [8] (for the comparison, the same 
bias and as is used). Our model is 2 times higher mainly 
due to our higher prediction for the IR galaxy emissivity. 
Similar results are found for other wavelengths. 



3. THE SIMULATIONS 

The simulations were computed by an IDL program that 
calculates the dark-matter power spectrum and spreads 
the galaxies in the map according to their correlation with 
the dark-matter density fleld. The CIB power spectrum is 



calculated as explained in Sect. 2. 

To create the maps, two assumptions were made: flrst 
that all the galaxies share the same spatial distribution 
independently of their luminosities; second that both IR 
and normal galaxies share the same spatial distribution. 
This second assumption was made to avoid too many 
free parameters in the simulations, the contributions of 
both populations being well separated in redshift this 
assumption is a weak one. 

The process for the creation of a virtual catalog can be 
summarised as follows. For a given wavelength, we create 
the map as a superposition of maps at different redshifts 
from z = to z — 6. The separation in redshift sHces 
decorrelates the emission from very distant regions of the 
modelled volume of the universe. In order to do so, we di- 
vided the maps in slices covering dz = 0.1. We can see the 
size of these sHces for different redshifts in Table [H For 
all redshift ranges the size of the slices is bigger than the 
measured comoving correlation lengths (for all populations 
of galaxies) . We then construct a brightness map for each 
redshift slice by adding: 1) a constant map with the mean 
surface brightness predicted by the LDP model for that z 
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Figure 8. Our power spectrum at 8 50 jim (continuous line) 
and that of iHaiman fc Knod l|2000D (dashed line). The dif- 
ferences between both models arise from the differences in 
the emissivities (see Fig.[T]). The lower level of the HK power 
spectrum comes from their lower emissivities. The emissivi- 
ties of HK are at lower z and therefore favour larger angular 
scales for the power spectrum relative to our model. 

Table 1. Physical size of the redshift slice dz = 0.1 (in 
Mpc) for different z. 



z 


1.0-1.1 


2.0-2.1 


3.0-3.1 


4.0-4.1 


Rdz=o.i (Mpc) 


233 


139 


93 


67 




Figure 9. Planck maps at 550 in MJy/sr with b=0 
(left) and b=1.5 {right). The maps simulate a region of the 
sky of 49 square degrees with 1024 pixels of 25 arcsec. 




Figure 10. Herschel maps at 500 /im in MJy/sr with b=0 
(left) and b=1.5 (right). The maps simulate a region of the 
sky of 0.3 square degrees with 1024 pixels of 2 arcsec. The 
small size of the maps makes it difficult to appreciate the 
effect of the large-scale clustering. 



Table 2. FWHM of the PSF for different wavelengths of 
observation (in arc seconds) for all the simulated maps. 

Wavelengths jfim) 350 550 850 1380 2097 
Planck HFI FWHM (") 300 300 300 330 480 ~ 

Wavelengths (fim) 250 350 550 " 

Herschel SPIRE FWHM (") 17 24 35~ 



were used for the simulations (5 — 0,0.75,1.5). Examples 
of maps at 500 fim (Herschel) and 550 //to (Planck) made 
with b=1.5 and b=0 are shown in Figs. [9] and [TOl The 
difference in the spatial correlations is easily noticed in the 
Planck simulations. On the other hand, the smaller size of 
the Herschel simulations makes it more difficult to see the 
correlation. 



slice, 2) a map of the fluctuations for the given bias pre- 
dicted by our spatial distribution model for that z slice. 
The fluctuations are not correlated between z sHces. The 
brightness map is then converted into flux map. At each 
luminosity, this can be converted into maps of numbers of 
sources. These numbers of sources are then redistributed 
into smaller z sHces (inside the 0.1 sHce) to refine the lu- 
minosity/flux relation. Note that all sources have the same 
underlying low frequency spatial distribution (but not the 
same positions) per dz = 0.1 slice. The position, luminos- 
ity, type (normal or starburst) and redshift of all sources 
are stored in a catalog. Since we know these four parame- 
ters for all the sources, we can now create maps of the sky 
at any given wavelength. To simulate the observations, the 
map is convolved with the point spread function (PSF) of 
the chosen instrument. 

For the purpose of this paper we have created 
Planck/HFI maps at 350, 550, 850, 1380 and 2097 microns 
and Herschel/SPIRE maps at 250, 350 and 500 microns. 
A description of the wavelengths and spatial resolution 
of the maps are given in Table [21 Three different biases 



The simulated maps and their asso- 
ciated ca talogs are publicly available at 
http:/ / www.ias.u-psud.fr/irgalaxies / simulations. php[ 

4. NOISE AND SOURCE DETECTION 

The simulations can be used to test the detection capa- 
bilities of Planck/HFI and Herschel/SPIRE. For the first 
time these simulations use an empirical model that repro- 
duces all the observational constraints from 5 fim to 1.3 mm 
and include the spatial correlation between the IR galax- 
ies and the dark matter density field for galaxies up to 
very low luminosities (L > 10^ Lq). They provide a useful 
tool for preparing future observations with Planck/HFI and 
Herschel/SPIRE. 

4.1. Detection of bright sources 

As stated previously bright sources dominate the power 
spectrum of the FIRB (see Fig. \3^. We therefore need 
to subtract them before studying the fiuctuations in the 
background. In this section we concentrate on detecting 
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them in three steps: 1) wavelet filtering, 2) detection, 3) 
measurement of the flux. 



— Wavelet filtering: Before trying to detect the sources 
we perform a wavelet transform of our simulated 
maps with the "atrou" algorithm. We remove spatial 
frequencies that are both higher and lower than the 
FWHM of the PSF. 

The small-scale filtering improves the estimation of 
the position of the sources when the instrumental 
noise is included in the simulations. In contrast to the 
confusion noise, the instrumental noise is not correlated 
for neighbouring pixels. This dominates errors in 
estimating the position of the sources. 

The large-scale filtering corrects for a bias in the 
detection algorithm. The algorithm searches for sources 
using the absolute value of the pixel and not its value 
relative to its environment. This biases the detections 
towards sources in bright regions. The removal of the 
large spatial fluctuations corrects this effect. 

The selection of spatial frequencies to be used for the 
detection has been manually optimised for each map 
to achieve a maximum number of reliable detections. 
This treatment is si milar to what w as done in the 
MIPS Spitzer maps l|Dole et all . l20Q4l l. A comprehen- 
sive study of the application of the wavelet filtering 
technique for the source detections at long wavelengths 
for Planck and Herschel/SPIRE is beyond the scope 
of this paper and has b een fully discussed in e.g. 
iLopez-Caniego et al.l 12003). 

— Detection algorithm: The algorithm is based on the 
"find" routine of the DAOPHOT library. In the filtered 
image the algorithm searches for peaks higher than a 
certain threshold Uthres- It uses the PSF shape and the 
neighbouring pixels to analyse whether the peak is the 
centre of a source. 

— Flux measurement: We developed a PSF fitting al- 
gorithm that we used in the original map (without 
filtering) to measure the fiux of the sources. We decided 
whether the detections are real or false by two criteria: 
1) proximity and 2) accuracy (see Sect. 4.1.1). 



4.1.1. Bad detections 

A detection is considered good or bad based on two 
criteria: 1) proximity with the position of an input source 
and 2) accuracy of fiux for this source. The former requires 
that our detection is closer than FWHM/5 to at least one 
"neighbour" source in our catalog. The latter requires that 
the difference between the fiux of one of the "neighbour" 
sources in the catalog and that of the detected source has 
to be smaller than the confusion and/or instrumental noise 
(see Tables [4] and [5]) . We consider the detection to be good 
only if both criteria are satisfied. 

The detection process also produces detections that do 
not comply with these criteria. We can see in Fig. [12] how 



the different (Jthres modify the rate of good-to-bad detec- 
tions. For a low detection threshold {(Jthres='2.<ymap, i-e. for 
example 290 mJy/pix at Planck 350 ^ra for a map with 
b=0 and no instrumental noise), the number of bad detec- 
tions can become bigger than that of real detections. For a 
higher detection threshold {athres=3o'map i-e. 440 mJy/pix 
at 350 /xm), we find that the good detections dominate the 
bad ones, but we do not detect as many faint sources. Thus 
the number of false detections depends strongly on (Jthres- 
For different scientific goals, it can be interesting to use dif- 
ferent (Jthres- For example, if we are interested in searching 
for objects at high redshifts, we could allow our detections 
to have 25% bad sources to be able to detect some interest- 
ing sources at high z. For studies of statistical properties of 
the sources, it would be necessary to use a stronger thresh- 
old. For our purpose, we used (Jthres = ScTmap (~ 10% of 

false detections). 

4.2. Instrumental and confusion noises 

Instrumental and confusion noises have been studied both 
separately and in combination in order to quantify their 
relative contribution to the total noise. The estimated 
instrumental noises per beam for Planck and Herschel are 
given in Table [3l The instrumental noise per beam for 
Planck is the average one over the sky for a 1-year mission. 
The instrumental noise per beam for Herschel is typical of 
large surveys. We take the sensitivity of the so-called level 
5 and level 6 of the Science Activity Group 1 (SAG 1) of 
the SPIRE guaranteed time team. 

We studied the standard deviation of the measured 
fiuxes in random positions for different maps. These 
maps were one of instrumental noise, three with dif- 
ferent bias (b=0, 0.75, 1.5) but without instrumental 
noise, and complete maps created by adding the map 
of instrumental noise to the three source maps. We call 
these maps hereafter instrumental-only, confusion-only, 
and complete-maps. We fit a Gaussian to the histogram 
of the fiuxes measured in these random positions and 
considered the standard deviation of this Gaussian as the 
best estimate of the standard deviation of the photometry 
of a source and therefore of the la instrumental, con- 
fusion, and total noise. Results are shown in Tables [4] and [H 

The confusion noise increases with the bias. This effect 
is noticeable for the Planck observations, but not for the 
Herschel ones because of the higher Herschel/SPIRE an- 
gular resolution. Also, for the considered Herschel/SPIRE 
surveys, the instrumental noise is always greater than the 
confusion noise. For Planck the correlation effect is more 
noticeable for longer wavelengths since they probe progres- 
sively higher redshifts and therefore higher dark-matter 
power spectra, as discussed in Sect 2.3 (see Fig. [5]). 

The total noise (Jc+i is close to the value 
cJ^+i = + '^1 (s^^ Tables [4] and [5|). For Planck at 
short wavelengths (350 fiui and 550 nm), the confusion 
noise is the dominant source of noise. The instrumental 
noise becomes dominant at 850/xm for b=0 and b=0.75. 
For longer wavelengths, it dominates for any bias. For 
Herschel the instrumental noise dominates the total noise 
for both the shallow and deep surveys. The confusion noise 
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Table 3. Simulation input instrumental noise per pixel of 
size equal to beam for Planck and for Herschel for a deep 
and a shallow survey . 



350 /<m 



Table 4. Noise on the retrieved sources with only instru- 
mental noise (cr/), confusion noise {(tc), and total noise 
{(Tc+i) in mJy for Planck/HFI. 



Wavelengths HFI (/im) 


350 


550 


850 


1380 


2097 




61.3 


39.2 


27.8 


16.7 


12.5 


ac b=0 


111.5 


41.5 


14.7 


4.6 


2.1 


ac b=0.75 


124 


54.3 


21.3 


7.7 


3.5 


ac b=1.5 


158 


79.8 


30.4 


10.9 


5.4 


ac+i b=0 


126.7 


62.3 


33 


17.4 


13.2 


(Tc+i b=0.75 


153.6 


75.3 


38.4 


19 


13.8 


ac+i b=1.5 


188.2 


95.3 


46.7 


21.1 


14.4 


Table 5. Instrumental noise (cr/), confusion noise (<Jc) and 


total noise {ac+i) in mJy for Herschel/SPIRE. 




Wavelengths SPIRE (/im 


250 


350 


500 






Deep tJi 


8.7 


11.3 


10.1 






Shallow ai 


15 


20.1 


18.2 






ac b=0, 0.75, 1.5 


4.6 


6.5 


5.5 






Deep ac+i 


9.8 


12.3 


11 






Shallow ac+i 


16 


20.8 


19 







is not strongly affected by the bias because of the small 
FWHM of the PSF. 



4.3. Completeness 

The first study of the point-sources detection limit for 
Planck was carried out bas ed on a generalisa t ion o f 
the Wie ner filter ing method ("Bouch et fc Gispertl . Il999f l. 
Recently iLopez-C aniego et al. ( 2006) used the most recent 
available templates of the microwave sky and extragalactic 
point source simulations, including both the radio and 
IR galaxies, to estimate the Planck detection limits. 
Here we revisit those results with new models for IR 
galaxies and the last noise estimates for Planck/HFI and 
Herschel/SPIRE. 

For the study of the completeness, a number Na of 
sources of equal fiux are randomnly distributed in the 
maps. Each source is placed far enough from the other to 
avoid these additional sources contributing to the confusion 
noise. The detection and photometry of these sources is 
carried out as described at the beginning of the section. We 
call Nq the number of good detections that comply with 
the "proximity" and "accuracy" criteria. The completeness 
for this flux Cf is then calculated as Cf — ^ x 100. 
The completeness of the detections of sources for a given 
flux depends on both the instrumental noise and the 
confusion noises. The results for the completeness are 
averaged over ^--^3000 individual fake sources per flux. An 



100 



Wavelengths HFI (^J.m) 


350 


550 


850 1380 2097 


cr/„st (mJy) 


31.30 


20.06 


14.07 8.43 6.38 


Wavelengths SPIRE (/im 


250 


350 


500 


ai„st Deep (mJy) 


4.5 


6.1 


5.3 


ainat Shallow (mJy) 


7.8 


10.5 


9.2 




300 



400 



500 
Flux mJy 



600 



700 



Figure 11. Study of completeness for Planck/HFI with 
6 = at 350/im. The horizontal straight line marks 
of completeness. 



350 fxm 



350 jLirr 




600 




Figure 12. Good detections (thick line) vs bad detections 
(thin line) Left: Histogram of good and bad detections using 
small (Jthres (290 mJy). Right: Histogram of good and bad 
detections using higher (Jthres (440 mJy) in the same map. 
Both plots have been done using a Planck simulated map 
with 6 = 0. 

Table 6. Completeness limits (in mJy) for the Planck/HFI 
maps with instrumental noise (C/), confusion noise (Cc) 
and both (Cc+/). We consider 6=0, 0.75, and 1.5. 



Wavelengths HFI (/^m) 


350 


550 


850 


1380 


2097 


Ci = 80% 


236 


157 


108 


67 


50 


Cc = 80% b=0 


516 


174 


60.5 


20 


8.6 


Cc = 80% b=0.75 


550 


239.5 


88.5 


30.5 


15.5 


Cc = 80% b=1.5 


684 


300 


121 


40 


24 


Cc+i = 80% b=0 


560 


234 


126 


71 


52 


Cc+i = 80%b=0.75 


607 


290 


141 


74 


55 


Cc+i = 80% b=1.5 


709 


360 


171 


80 


58.5 



example of the completness at SSO/xto is shown in Fig. 
[TTl Results for all wavelengths are given in Tables [6] and 
[7l They are consistent with the instrumental, confusion, 
and total noises given in Sect. 4.2 and the conclusions 
from that section remain valid for the completeness. For 
simulated maps including both extragalactic sources and 
instrumental noise, we find that the 80% completeness 
level coincide with fiux limits around 4 — 5(7. 



Taking these 80% completness limits as a detection 
threshold, the prediciton for the number of sources de- 
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Table 7. Completeness limits (in mJy) for the 
Herscliel/SPIRE maps with instrumental noise (C/), con- 
fusion noise (Cc), and both (Cc+/). 



Wavelengths SPIRE(^im) 


250 


350 


500 


Deep Cj = 8 


0% 


33 


45.9 


37.9 


Shallow Ci = 


80% 


57.3 


84.2 


67.3 


Cc = 80% 


35 


32 


27.4 


Deep Cc+i = 


80% 


49.8 


64.4 


48.4 


Shallow Cc+i - 


- 80% 


70 


96.5 


75.5 



tected directly by Herschel/SPIRE for the deeper survey 
considered here is 8.3x10 /sr at 250 ^m, l.lxlO^/sr 
at 350 ^m, and 1.8xl0^/sr at 500 /im. The fraction of 
resolved CFIRB varies between 8 and 0.3% from 250 to 
500 /im. 



5. CIB fluctuations 

The Planck and Herschel/SPIRE surveys allow an unprece- 
dented search for CFIRB fluctuations associated with large- 
scale structure and galaxy clustering. Background fluctua- 
tions probe the physics of galaxy clustering over an en- 
semble of sources, with the bulk of the signal contribution 
originating from sources well below the detection thresh- 
old. Thus a comprehensive fluctuation analysis is an es- 
sential complement to the study of individually detected 
galaxies. In this section, we restrict ourselves to predictions 
for Planck/HFI, excluding the 143 and 100 GHz channels. 
At these low frequencies, we are dominating by the non- 
thermal emission of the radiosources - that are not includ- 
ing in our model - and the Poisson terrn dom inates the 
clustering term (e.g. lGonzalez-Nuevo et aL 1. 1200,1 ). Also, we 
exclude the Herschel/SPIRE case since our simulations in- 
clude the clustering of CIB sources in two different halos 
(2/i), but not the clustering within the same halo (l/i). The 
Ih term dominates for t > 3000 and will be accurately 
measured by Herschel/SPIRE. Only large-scale surveys can 
put strong constraints on the 2h term. Measuring the 2h 
clustering with CFIRB anisotropics is one of the goals of 
Planck/HFI. 

5.1. Contributors to the angular power spectrum 

From the far-IR to the millimeter, the sky is made up of the 
CFIRB and two other sources of signal, the galactic cirrus 
and the CMB (we neglect the SZ signal). Understanding 
our observations of the CFIRB requires understanding the 
contributions from these two components which act for us 
as foreground and background contamination. 

The galactic cirrus acts as foreground noise for the 
CFIRB. The non- white and non-Gaussian statistical prop- 
erties of its emission make it a very complex foreground 
component. The power spectrum of the I RAS 100 /xm 
emis sion i s characterised by a power law iGautier et al.1 
(e.g. Il993 l. Here we compute the angular power spectrum 
of th e dust emission following iMiville-Deschenes et al.l 
l|2007( ). These authors analysed the statistical properties 
of the cirrus emission at 100 fiui using the IRAS/IRIS 
data. We used their power spectrum normalization and 
slope (varying with the mean dust intensity at lOO/xm). 



Using the average \b\ > 30° spectrum of the Hl-correlated 
dust emission measured using FIRAS data, we converted 
the 100 /zm power spectra to the Planck wavelengths. For 
the discussion, we considered both the total cirrus fluctu- 
ations and 10% residual fluctuations. These 10% could be 
achievable with Planck in low dust-column-density regions 
containing ancillary HI data. 

The CMB acts as background noise for the CFIRB. 
However, the CMB angular power spectrum is known 
to an accuracy of 1% or better (Planck-HFI web site: 
http://www.planck.fr/). This combines with its well- 
known spectral dependence to allow for a clean subtraction 
of its contribution that in turn allows for detections of 
the CFIRB C/ even for wavelengths where the CMB 
dominates. We consider for the rest of the discussion a 
conservative assumption, that is that the residual CMB 
fluctuations are approximately 2%. 

The angular power spectrum of IR galaxies is composed 
of a correlated and a Poissonian part. As discussed in Sect. 
2.1, the contribution to the Poissonian part is dominated 
by relatively faint sources after subtracting the brightest 
galaxies (see Fig. [3]). We consider that we can remove 
sources brighter than our 80% completeness detection Hmit 
(see Tables [6] and III) . For doing so, we use the technique 
described in Sect. 4.1 for measuring the position and fluxes 
of the sources and once these are known we subtract a 
PSF with the measured flux from the map. 

The correlated Ci is obtained as described in Sect 2.1. 
For the rela tive error on the power spectrum, we follow 
iKnoxl ||1995( ): 

Ci \aJ \2l + lJ y ^ NCiWi J 

where A is the observed X the rms noise per 

pixel (instrumental plus confusion), N the number of 
pixels, and Wi the window function for a map made with 
a Gaussian beam Wi = e~' . 



5.2. Detectability of CFIRB correlated anisotropies 

The study of the Ci on difi^erent scales allows us to study 
different aspects of the physics o f the environment of 
IR galaxies (see ICoorav fc Shethl . l2002f ). Large scales 
{I < 100) give information on the cosmological evolution 
of primordial density fluctuations in the linear phase and 
therefore on the cosmological parameters. Intermediate 
scales are mostly influenced by the mass of dark halos 
hosting sources, which determines the bias parameter. 
Small scales {I > 3000) probe the distribution of sources 
in the dark-matter halos and therefore the non-linear 
evolution of the structures. This non linear evolution was 
not accounted for in our model. 

We can see in Fig. [13] the different contributions to the 
Ci for Planck/HFI: the correlated CFIRB C; with 6=1 
and 6 = 3 with their respective error bars {A£/£=0.5), 
Poissonian fluctuations, dust, and CMB contributions. 
The dust and CMB Ci are plotted both before and after 




Figure 13. Power spectra for Planck of the correlated IR galaxies with b=l (thick continuous Hne) and b=3 (thin 
continuous Hne), Poisson fluctuations for sources fainter than 709, 363, 171, 80 mJy at 350, 550, 850 and 1380 fim 
respectively (short dashed line), dust with an HI column density of 1.5 10^" at/cm? and 2.7 10^° atjcrv? for the 400 
Sq. Deg. and ~ 15000 Sq. Deg. considered regions respectively (dotted line), dust/ 10 (dotted-dashed Hne), CMB (dotted 
Hne) and CMB/50 (dashed line). Error bars are computed using A^/£=0.5. 



Table 8. Multipole I ranges where the CFIRB correlated anisotropies are higher than the cirrus and CMB components. 



Wavelengths (/im) 


350 


550 


850 


1380 


400 Sq. Deg. b=l 
400 Sq. Deg. b=3 
15000 Sq. Deg. b=l 
15000 Sq. Deg. b=3 


600 < ^ < 1200 
120 < £ < 3500 

Indetectable 
280 < ^ < 3500 


200 < ^ < 1800 
75 < £ < 5000 
550 < ^ < 1800 
130 < ^ < 5000 


1000 < 
240 < L 
1200 < 
230 < L 


I < 2200 

< 6000 
£ < 2200 

< 6000 


Indetectable 
1800 <£< 6500 

Indetectable 
1800 <£< 6500 
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Table 9. Multipole I ranges where the CFIRB correlated anisotropies are higher than the residual cirrus and CMB 
components (10%, and 2% respectively). 



Wavelengths (/xm) 


350 


550 


850 


1380 


400 Sq. Deg. b=l 


200 < Z < 1200 


80 < / < 2100 


60 < ; < 2200 


Indetectable 


400 Sq. Deg. b=3 


40 < Z < 3500 


30 < / < 5000 


25 < ; < 6000 


m<l< 6500 


15000 Sq. Deg. b=l 


250 <; < 1200 


130 < 1800 


100 < Z < 2200 


Indetectable 


15000 Sq. Deg. b=3 


80 < Z < 3500 


55 < / < 5000 


45 < Z < 6000 


80 < Z < 6500 



implementing the corrections discussed above (10% and 
2% residuals, respectively). For the dust we selected for 
the estimation of the column density a field centred on the 
SWIRE/ELAIS SI field. This field has a very low level 
of dust contamination. For a 400 sq. deg. area centred 
at {l,h) — (311°,— 73°), we have an average HI column 
density of only 1.5 10^*^ at/cm^. We also take the whole 
sky above \h\ > 40° (~ 150000). The average HI column 
density is 2.7 10^" at/cm^. 



At all wavelengths, the correlated CFIRB C; is dom- 
inated by the other components on both very large and 
very small scales. The angular frequencies where the C; 
is dominated by the galaxy correlation depends mainly 
on the wavelength, if the bias and cirrus level are fixed. 
On large scales the C^^^^^ is dominated by the dust for 
wavelengths up to 550 ^m, and for longer wavelengths it 
is dominated by the CMB. On small scales the Poissonian 
fiuctuations dominate the power spectra. 



We can see in Tables [8] and [9] the ranges for which 
CFIRB-correlated anisotropies are dominating depending 
on whether we consider a partial subtraction of the cirrus 
and CMB or not. On small angular scales (large £), the 
range of € where the correlated Ci dominates increases 
with the wavelength. On large angular scales (small £), the 
results depend on whether we consider partial subtraction 
of the cirrus and CMB or not. 



Table [8] shows the I ranges for where CFIRB correlated 
anisotropies are dominating when no subtraction of the 
dust or CMB has been performed. At large angular scales 
and wavelengths up to 850 ^im the detectability is better 
for the small map because of its lower N(HI) column 
density. Without any dust contamination correction, this 
prevents us from taking advantage from the smaller errror 
on the estimation of the Ci in the very large area surveys 
and makes both kind of observation complementary. In 
clean regions of the sky, the best wavelengths for the 
observation of the large scales CFIRB anisotropies are 350 
lim and 550 /zm. 



The large-scale detectability drastically improves in the 
case of partial subtraction of the CMB and dust cirrus as 
seen in Table [9l As discussed above we expect to be able to 
subtract the dust fiuctuations and the CMB to a 10% and 
2% of their original level. This does not change our ability 
to measure the correlated CFIRB C; on small scales (due to 
the high confusion noise), but will allow us to probe larger 
scales. 



6. Summary 

This paper presented new simulations of the cosmic in- 
frared and submiUimeter background. The simulations are 
based on an empirical model of IR galaxy evolution (the 
LDP model) combined with a simple description of the 
correlation. The IR galaxy spatial distribution follows that 
of the dark-matter halo density field, with a bias parame- 
ter accounting for the possibility of the luminous matter 
being more or less correlated with the dark matter. The 
simulated maps and their catalogs are publicly available 
at http:/ /www. ias.u-psud.fr/irgalaxies/simulations.phpl 
Other maps are available upon request. These maps are 
intended to be a useful tool for planning future large IR 
and submiUimeter surveys. In this paper, we used the maps 
to predict the confusion noise and completeness levels 
for Planck/HFI and Herschel/SPIRE. We also predicted 
the power spectra of correlated CFIRB anisotropies for 
Planck/HFI which will be a major advance in the study of 
the CFIRB anisotropies at large scales (i.e. I < 2000 - 5000 
depending on the wavelengths). Further analysis of the 
CFIRB anisotropies, including the use of stacking analysis 
to isolate the anisotropies in different redshift ranges, will 
be presented in a second paper, now in preparation. 
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